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Abstract 

We study disordered spin-1 quantum chains with random exchange 
and biquadratic interactions using a real space renormahzation group 
approach. We find that the dimerized phase of the pure biquadratic 
model is unstable and gives rise to a random singlet phase in the 
presence of weak disorder. In the Haldane region of the phase diagram 
we obtain a quite different behavior. 

KEYWORDS: A. Disordered Systems; D. Phase Transitions, Thermody- 
namic Properties. 



Although much progress has been made on the understanding of the phys- 
ical properties of ordered isotropic quantum spin chains, for both integer and 
half-integer spins ^ , much less has been achieved in the random case in par- 
ticular for spin-1 chains. For the spin- 1/2 random exchange Heisenberg anti- 
ferromagnetic chain (REHAC) theoretical advance ^ has been motivated 
by experimental results mainly in organic materials IQ. The theory in this 
case gives a good description of the physical properties of these systems which 
are governed by low energy excitations [^. Also recently an extensive study 
using a renormalization group approach has led to a deeper understanding of 
the zero temperature phase diagram of general spin-1/2 disordered quantum 
chains On the other hand the case of random, isotropic, spin-1 chains 
has only recently been addressed both experimentally 0, |^ and theoretically 



In this Communication we study a one-dimensional, spin-1, Heisenberg 
system with random exchange and biquadratic interactions. We use a real 
space renormalization group approach which has proved to be very sucessful 
for spin-1/2 quantum chains 0, ^. We discuss the ground state properties 
and extend the approach to finite temperatures to obtain the thermodynamic 
properties. 

We consider the general bilinear, biquadratic, spin-1 chain with L spins 
described by the Hamiltonian 
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where Jr and Dr are random nearest-neighbor interactions with probabihty 
distributions Pj{Jr) and PoiDr) respectively and such that < Jr < J and 
< Dr < D. The Sr are spin-1 operators and we are interested in the hmit 
L ^ oo. 

We study in this Communication the effect of disorder in two distinct 
situations: 

1. The case where the biquadratic is the dominant interaction such that 
the cutoff D >> J. We recall that the ordered system with Dj. = D 
and furthermore J,. = corresponds to the so called KBB point of the 
general phase diagram of the Hamiltonian given by Eq.l [|]. At this 
point the purely biquadratic Hamiltonian can be solved exactly yielding 
a dimerized phase characterized by a two-fold degenerate ground state 
with a finite gap to the excited states. 

2. The cutoff J of the exchange interaction is much larger than that of the 
biquadratic couplings, i.e. J » D. This corresponds to investigate 
the effect of disorder on the section of the phase diagram of the pure 
system where a Haldane gap occurs . The ground state of the ordered 
system {Jr = J, Dr = 0) in this case is unique and also has a gap for 
excitations. 

In order to deal with the Hamiltonian, Eq.l, we generalize, for the case 
of spin-1 and biquadratic interactions, a real space renormalization group 
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method introduced by Ma, Dasgupta and Hu to treat the spin-1/2 RE- 
HAC. The usefulness and hmitations of this method have been intensively 
discussed by Fisher in recent papers on quantum random spin-1/2 chains. 
In at least one of the situations we study here this approach converges to the 
exact solution of the random problem. The method consists of identifying 
the pair of spins with the strongest coupling in the random chain and elim- 
inating it by considering the interaction with the neighboring spins of this 
pair as a perturbation. This procedure generates new couplings J' and D' 
between the two spins neighboring the eliminated pair and consequently the 
form of the distribution Pj{Jr) and Po{Dr) are modified. One also gains a 
contribution to the free energy of the system from the eliminated pair. This 
process is then successively applied until a single pair of spins remains. 

We start at zero temperature and consider the case of purely biquadratic 
interactions and no exchange. In the ordered mentionned above this 

corresponds to the so-called KBB point |^ of the general phase diagram of 
the spin-1 chain given by the Hamiltonian of Eq.l. We assume a distribution 
of random biquadratic couplings and initially take J,. = 0. Generalizing the 
procedure of Ma et al. for this problem, we identify the pair of spins with 
the strongest biquadratic coupling. When we apply the elimination transfor- 
mation to this strong coupled pair we find an effective coupling between the 
spins neighboring the eliminated pair which is is given by: 

D' = (2) 

9 D ^ ' 



where Di and D2 are the bonds connecting the strongest coupled pair of 
spins, with interaction strength D, to their neighbors. It is clear that in this 
case the successive application of the elimination transformation generates 
weaker and weaker couplings and the distribution of biquadratic interactions 
becomes peaked close to the origin as the cutoff D decreases. The physical 
properties of the system are then dominated by low energy excitations and 
this new phase is very similar to the random singlet phase of the disordered 
spin-1/2 random antiferromagnetic chain 0. Consequently disorder has a 
dramatic effect at the KBB point since in the pure case the system has a 
dimerized phase with a gap for excitations and a short correlation length. The 
introduction of weak randomness on the biquadratic couplings gives rise to 
a new disordered phase with low energy excitations and long-range magnetic 
correlations. This behavior is universal in the sense that independently of 
the form of the initial distribution of couplings Pd(D^.) this converges to a 
fixed point distribution, which is approximately described by a power law 
with a singularity at the origin 0. Notice that the equation we find here 
renormalizes to weak coupling faster than in case of the spin-1/2 chain where 
the pre-factor of the recursion relation equivalent to Eq.2 is 1/2 p|. We have 
found here a rather unusual situation where a phase with a gap and a short 
correlation length has become completely unstable and furthermore with 
long-range magnetic correlations due to the introduction of arbitrarily weak 
disorder p. 
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We include now a weak random exchange interaction between nearest 
neighbor spins such that the cutoff J « D. A generahzation of the ehmi- 
nation transformation gives, for the new exchange and biquadratic couphngs 
between the spins neighboring the ehminated pair, the following coupled re- 
lations 

^, 2 D1D2 

D' = (3) 

9 D + Jo ^ ' 

J' = i -^'-^^ 



3 {3D + Jo) 

where as before Di, D2, J\ and J2 are the bonds connecting the strongly 
coupled pair of spins, with interactions D and Jq between them, to their 
neighbors. Since D » J we expect, as before, that the renormalized distri- 
butions converge rapidly to a weak coupling situation, i.e., to develop more 
and more weight close to the origin as the cutoff decreases. Again the system 
attains a random singlet phase dominated by low energy excitations. 

We remark that the recursion relations for the biquadratic interactions 
iterate to zero faster than that of the exchange couplings. It may happen 
that, depending on the starting distributions, the physics at sufficiently low 
energies will be dominated by the exchange interactions. In this case it is 
convenient to define the strongest coupled pair as that for which 3D -|- J, the 
energy to create an excited state for a pair, is larger. 

For completeness we give now the expression for the energy associated 
with the eliminated pair and from which we can compute the total ground 
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state energy 



E' ^ -2Jo - 4D - -{D, + D2) 



3 Jo + 3D 



8 Df + Dl- 



(4) 



27 Jq + D 



We conclude this analysis of case 1, i.e. D >> J, pointing out that 
we have obtained a new disordered phase for the random biquadratic spin-1 
Hamiltonian which is similar to the random singlet phase of the spin- 1/2 



Case 2) We start with the situation where Dr — 0. For the random 
exchange spin-1 chain, the elimination procedure yields the following relation 
for the new coupling between the spins neighboring the eliminated pair 



We find here quite a distinct situation from the previous case, since this 
recursion relation appears to iterate to strong coupling. The generation of 
couplings which are larger than those eliminated would in fact invalidate the 
present approach based on perturbation theory. Contrary to the previous 
cases, now the results depend strongly on the form of the starting distribu- 
tions. When we include the biquadratic interaction, ior D << J, the new 
recursion relations are similar to those obtained before, namely 



REHAC. 



j' = ; 



4J1J2 



(5) 



3 J 




4 J1J2 
3(J + 3L>o) 



(6) 
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^> _ 2 D1D2 
9J + D0 

ioT D << J. It is clear that the biquadratic is an irrelevant interaction, in 
the renormalization group sense, for the random exchange spin-1 chain since 
the starting distribution P£,{Dr) iterates to a distribution with most of it's 
weight at small coupling much faster than the exchange distribution. We 
notice that the biquadratic interactions reduce the new exchange couplings 
generated by the elimination transformation. Whether this is sufficient to 
guarantee the validity of the present approach will depend on the nature of 
the starting distributions and the values of the cutoffs. 

In order to get a better understanding of the random exchange spin-1 
chain, we have generalized the elimination transformation procedure to finite 
temperatures |]^. In the case = we get after some lenghty calculations 

J' = l^wm (7) 

where Ji and J2 as before are the bonds connecting the strongly coupled pair 
of spins to their neighbors. We point out that this transformation does not 
generate new terms preserving the form of the original Heisenberg Hamil- 
tonian. The free energy has a contribution from the eliminated pair which 
turns out to be, 

F' = F,- l-^^Vm (8) 

with 

Fo = -2J - \Ln{l + Se"^'^ + Se"^^-^)) (9) 
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The quantities W(y) and V(y) are given by |0: 

.... . l-(3/8)e-ni+y)- (5/8)6-3^1 + 3y) 

^^y^ = 1 + 3e-y + 5e-3. (1°) 

and 

... ^ l-(3/8)e-^(l-y)-(5/8)e-^^(l-3y) 

^^y^ = 1 + 3e-y + 5e-3. ^^^^ 

Notice that for ksT > 0.37J the function W{PJ) < (3/4) and this guaran- 
tees the vahdity of the present approach for sufficiently high temperatures 
independent of the starting distribution. 

For certain probabihty distributions, which we associate with the case 
of strong disorder , the successive ehmination transformations give rise to 
weaker and weaker couphngs as the cutoff J decreases. In this case, for 
sufficiently small J, the distribution Pj{Jr) becomes peaked at ~ and 
can be approximated by a power law ||, |^ , 

Pj{x, J) ^ -x-'+" (12) 

where x = Jr/J, J is the actual cutoff and the factor {a/ J) arises from 
the normalization condition. This distribution leads to a low temperature 
behavior of the thermodynamic quantities which can be described by power 
laws 0,0], 

F oc = _t|^ oc T-^^ (13) 



and 



XocT~^+^%mocif^'' (14) 



where 7c is the specific heat exponent, 7^ and 7/1 are related to the exponent 
a of P(x,J) and of the free energy F . The susceptibihty exponent is given, to 
a good approximation by, 7^ ~ la and furthermore 7/1 ~ 7^ H and m are 
respectively the uniform external field and magnetization. Since the power 
law distribution is not the exact fixed point distribution for the elimination 
transformation, a and consequently 7c, 7s and depend on temperature 
and cutoff. For initial power law distributions, Eq.l2, this dependence is 
very weak [@, 0]. 

The results contained in Eqs.l3 have been explicitly verified for the strong- 
ly disordered cases of initial uniform, P{K) = 0(1 — K), and power law 
distributions (Eq.l2) The free energy was obtained numerically from 
the accumulation point of Eq.8 for successive elimination transformations. 
Then Eqs.l3 yield the exponent a and the specific heat 0. For the same 
starting distribution and temperature range we found that the exponent a 
for the spin-1 chain is consistently larger than that for the spin-1/2. In fact, 
for any initial distribution which attains a power law form whenever the 
cutoff is sufficiently small, we can show that the exponent a decreases more 
rapidly when the cutoff is further reduced for S = 1/2 than for 5 = 1 [^]. 
These theoretical results are in agreement with experiments in the strongly 
disordered, one-dimensional, spin-1 system MgYOBO-^ which shows power 
law behavior of the susceptibility in a wide temperature range although 
less singular than in the isostructural spin-1/2 material MgTiOBO^ 0. 
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The elimination transformation for S = 1, Eq.7 or Eq.5, applies for 
strongly disordered chains , which we generically define as those characterized 
by initial distributions which yield thermodynamic properties at low temper- 
atures as described by Eqs.l3 and 14. For these distributions couplings larger 
than the one eliminated are generated with negligeable probability and the 
cutoff decreases rapidly enough to yield meaningful thermodynamic behav- 
ior as in the cases of uniform and power law distributions. Otherwise our 
results represent a good approximation for the high temperature behavior of 
disordered spin-1 chains. 

We have extended the elimination transformation method for a quantum 
Heisenberg chain with arbitrary spin. We find the following recursion relation 



This result shows that perturbation theory fails for large spins since the 
excited states become nearly degenerate with the ground state. The critical 
value of S is between S — 1/2 and S — 1 so that a universal random singlet 
type of phase in the Heisenberg chain with exchange interactions only is 
guaranteed for spin- 1/2. In this case the quantum fiuctuations are sufficiently 
strong and the renormalization group equations yield a fiow to a fixed point 
universal distribution singular at the origin. 

In conclusion we have generalized an elimination transformation proce- 
dure to study the effect of disorder on quantum antiferromagnetic chains 
with exchange and biquadratic interactions. In case the biquadratic inter- 



J' = '^S{S + 1) 



J1J2 



(15) 



J 
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action dominates, we find a new disordered phase for spin-1 cliains wliicli is 
essentially a random singlet phase. We find that weak disorder destroys the 
dimerized phase of the ordered system giving rise to low energy excitations 
and long range correlations. In the Haldane region of the phase diagram of 
the pure chain our results indicate that weak disorder has not a profound 
effect although the present approach fails in this case at least at low tem- 
peratures. However for starting distribuitons which are strongly disordered 
we obtain a thermodynamic behavior which is similar to the random singlet 
phase of disordered spin- 1/2 REHAC but with weaker singularities. 
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